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Abstract:  Generalized  polynomial  chaos  (gPC)  and  the  probabilistic  collocation  method  (PCM) 
are  finding  considerable  application  to  problems  of  interest  to  engineers  in  which  random  parame¬ 
ters  are  an  essential  feature  of  the  mathematical  model.  So  far  the  applications  have  been  mainly 
to  stochastic  partial  differential  equations,  but  we  extend  the  method  to  volume-integral  equations, 
which  have  met  great  success  in  electromagnetic  nondestructive  evaluation  (NDE),  especially  with 
eddy-currents.  The  problems  of  main  interest  to  the  NDE  community  in  this  connection  are  con¬ 
cerned  with  the  issue  of  ’propagation  of  uncertainty’  when  the  relevant  parameters  are  not  well 
characterized,  or  are  known  simply  as  random  variables.  We  demonstrate  the  ideas  by  considering 
a  metallic  surface  that  has  undergone  a  shot-peening  treatment  to  reduce  residual  stresses,  and 
has,  therefore,  become  a  random  conductivity  field. 

Keywords:  volume-integral  equations,  electromagnetic  nondestructive  evaluation,  generalized 
polynomial  chaos,  probabilistic  collocation  method 

1.  Introduction 

Critical  structural  components  are  often  treated  by  shot-peening  or  low-plasticity  burnishing  to 
reduce  residual  tensile  stresses.  This  process  often  leaves  the  surface  of  the  component  rough  with 
a  profile  that  can  only  be  defined  as  a  random  process.  It  is  necessary  to  characterize  the  surface, 
and  the  most  common  method  is  through  the  use  of  eddy-current  nondestructive  evaluation  (NDE). 
This  method  measures  the  conductivity  variation,  which,  because  it  simulates  the  rough  surface, 
becomes  itself  a  random  field.  This  paper  describes  an  approach  to  modeling  electromagnetic 
interactions  with  a  random  conductivity  field  in  the  context  of  NDE.  The  computational  engine  is 
based  on  stochastic  volume- integral  equations. 
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2.  Volume-Integral  Equations 


We  have  introduced  the  notion  of  volume-integral  equations  and  their  application  in  eddy-current 
NDE  in  a  recent  series  of  papers  [1-3].  We  will  remind  the  reader  that  the  essential  feature,  for 
our  purpose,  is  that  the  ’anomalous  region’  (the  ’flaw’)  is  decomposed  into  a  regular  grid  of  voxels, 
each  of  which  has  uniform  conductivity.  Because  the  surface  that  this  grid  simulates  is  random,  so, 
too,  is  the  conductivity  variation  from  voxel  to  voxel.  In  other  words  we  have  defined  a  random 
conductivity  field,  with  each  voxel  conductivity  being  a  random  variable.  The  discretization  process, 
of  course,  transforms  an  infinite-dimensional  field  (the  ’true  surface’)  into  a  finite-dimensional  field, 
but  still  with  a  very  large  number  of  random  variables.  It  is  this  feature  that  presents  a  numerical 
challenge  in  modeling  stochastic  integral  equations. 

3.  Theoretical  Background 

If  there  are  N  voxels,  then  the  impedance  computed  by  the  volume-integral  equation  can  be  written 
Z(l,a i,---,  <7jv),  where  l  is  the  position  of  the  probe  coil,  and  the  average  value  of  the  impedance 
can  be  written 

/OO  POO 

■  ■  ■  ■  ,aN)p(ai,  ■  ■  ■  ,aN)dai  ■  ■  ■  daN  ,  (1) 

-oo  J —oo 

where  p(a\,  ■  ■  ■ ,  <jn)  is  the  joint  probability  density  of  the  conductivity  held.  N  will  generally  be 
large,  so  we  must  consider  a  mapping  into  a  lower-dimensional  random  space.  One  such  mapping 
is  the  Karhunen-Loeve  expansion  [4,5].  The  discrete  version  of  the  K-L  expansion  states  that  every 
random  vector,  V,  has  a  decomposition 


N 

V  =  K£nVn  ,  (2) 

n=  1 

where  {An}  and  {vn}  are  the  eigenpair  solutions  of  the  covariance  matrix  of  V,  and  {£n}  are 
uncorrelated  random  variables,  =  8mn.  It  is  hoped  that  this  expansion  converges  after  M 

terms,  where  M  «  N  in  (1).  The  {CnjnLi  then  become  the  new  random  variables  in  expressions 
such  as  (1). 

We  will  apply  the  K-L  expansion  to  the  double-exponential  measured  correlation  function, 
C(x,x')  =  exp(— \x  —  x'\/L),  which  arises  from  shot-peened  copper  for  two  different  peening  inten¬ 
sities  [6].  The  discrete  covariance  matrix,  Grnn ,  for  this  function  is  given  by  [7]: 


Gran  =  -2L2  (l  -  e~s/L  -  8 / l)  , 

m  =  n 

=  -2 L2e-(n-m),5/L  _  cosh(<5/L))  , 

m  <  n 

(3) 

=  -2 L2e-(m-nWL  (1  -  cosh (8/L)) , 

m  >  n  , 

where  L  is  the  correlation  length  of  the  process,  and  5  is  the  length  of  a  one-dimensional  cell. 

Figure  1  shows  the  normalized  eigenvalue  spectrum  for  a  32  x  32  matrix  for  three  values  of  the 
ratio,  L/8,  and  sample  functions  for  the  anomalous  conductivity  profiles  are  shown  in  Figures  2-4 
for  these  same  values  of  L/8.  Sample  functions  for  the  case,  L/8  =  0.1,  are  shown  in  Figure  2,  while 
Figure  3  shows  three  sample  functions  for  the  conductivity  profile  when  L/8  =  1,  and  Figure  4 
shows  three  sample  functions  for  the  conductivity  profile  when  L/8  =  3. 
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Normalized  eigenvalues  as  a  function  of  L/delta 


Figure  1:  Normalized  eigenvalue  spectrum  for  the  double-exponential  covariance  function  for  three 
values  of  the  critical  ratio,  L/5. 


Three  conductivity  profiles  for  L/delta=0.1 


Figure  2:  Three  sample  functions  for  the  conductivity  profile  when  L/5  =  0.1.  These  functions  are 
the  departure  from  the  mean  value  of  Ohost  =  3.02  x  105  S/m.  We  assume  a  uniform  probability 
density  function,  centered  at  zero  and  with  variance  =  1,  for  the  random  variables,  {£*},  in  the 
Karhunen-Loeve  expansion. 
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Three  conductivity  profiles  for  L/delta=1 


Figure  3:  Three  sample  functions  for  the  conductivity  profile  when  L/5  =  1.  These  functions  are 
the  departure  from  the  mean  value  of  <7host  =  3.02  x  105  S/m.  We  assume  a  uniform  probability 
density  function,  centered  at  zero  and  with  variance  =  1,  for  the  random  variables,  {£*},  in  the 
Karhunen-Loeve  expansion. 


Three  conductivity  profiles  for  L/delta=3 


Figure  4:  Three  sample  functions  for  the  conductivity  profile  when  L/5  =  3.  These  functions  are 
the  departure  from  the  mean  value  of  rJhost  =  3.02  x  105  S/m.  We  assume  a  uniform  probability 
density  function,  centered  at  zero  and  with  variance  =  1,  for  the  random  variables,  {£*},  in  the 
Karhunen-Loeve  expansion. 
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4.  gPC  and  PCM 


We  will  assume  that  the  new  random  variables,  {£1,  •  •  •  ,£m},  that  are  the  result  of  the  K-L  ex¬ 
pansion,  are  independent,  identically  distributed  random  variables  (iid),  with  the  common  density 
function,  7 r(£).  The  generalized  polynomial  chaos  (gPC)  expansion  of  degree  N  for  Z(l,£  1,  •  •  •  ,£m) 
is  given  by  [5]: 

■  ■  ■  >£m)  =  E  ,  (4) 

0<\i\<N 

where  |i|  =  i\  +  ■  ■  ■  +im  and  {</>fc(£j)}Eo  the  se^  °f  univariate  gPC  basis  functions  in  C  of  degree 
0  <  k  <  N.  They  are  orthogonal  polynomials  associated  with  the  density  function,  7r(£),  in  the 
sense  that 

=  J  =  LnTm.  0  <  to,  n  <  iV  ,  (5) 

where  7 m  is  a  normalizing  constant. 

With  this  expansion  in  hand,  we  can  calculate  various  statistical  properties  of  Z(l,£  1,  •  •  •  ,£m) 

[5].  For  example,  [Zq{1)]m  is  the  average  value  of  Z .  and  the  variance  and  covariance  are  given  by: 

Var[Z(Z,  6, ' ' ' ,  6f)]  =  E  7ii'"  7 iMZ\  (h)  •  •  •  Z?M  0 h ) 

0<|i|<iV 

Cz{h,h)  =  E  7u  ■  ■ ■  ■  ■  ZiM(li)Zh(l2)  ■  ■  ■  ZiM(l2)  ,  (6) 

0<|i|<iV 

where  the  sums  exclude  |z|  =  0. 

5.  HDMR  and  ANOVA 

High-dimensional  model  representation  (HDMR)  and  analysis  of  variance  (ANOVA)  seek  to  reduce 
the  complexity  of  problems  with  a  large  number  of  dimensions,  and  are  the  subject  of  considerable 
contemporary  research  in  computational  methods  [8-11]. 

Instead  of  the  gPC  expansion  of  (4),  we  consider  an  ANOVA  decomposition  as  (we  suppress 
the  scan  variable,  l,  to  simplify  notation): 

M  M 

Z(£i,  ■  ■  ■  ,€m)  =  Zq  +  E  Zjx (£ji)  +  E  Zj1j2(^j1,Cj2) 

h  ji<j2 

+  E  Zj  1J2  43(^1  >  C?2 >  C73)  T  •  •  •  +  Zjlt...tjM  (Oi  !  '  ■  ■  )  £,jM  )  ,  (7) 

jl<j2<j3 

where  Zo  is  a  constant  function,  the  { Zn  }  are  one-dimensional  functions,  {Zj1j2}  are  two-dimensional 
functions,  and  so  on,  yielding  2A1  different  terms  [9].  The  various  functions  are  orthogonal  to  each 
other  (stochastically  uncorrelated),  which  means  that  Zq  is  the  mean  value  of  Z,  and  the  remaining 
terms  are  partial  variances,  whose  sum  is  the  total  variance  of  Z. 

ANOVA  and  PCM  will  be  used  to  compute  the  stochastic  response  of  the  volume-integral 
equation  given  random  parameters  of  the  anomalous  region  that  evolves  from  shot-peening  or 
other  surface  treatment  methods.  This  response  will  be  expressed  in  terms  of  average  values 
and  variances,  or  even  probability  densities.  These  will  have  further  ramifications  in  computing 
probability-of-detection  (POD)  for  various  testing  setups,  as  well  as  inverting  measured  data  to 
estimate  flaw  parameters,  etc. 
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